Search for: All records

We present a robust methodology for identifying photometric binaries in star clusters. Using Gaia DR3, Pan-STARRS, and Two Micron All Sky Survey data, we self-consistently define the cluster parameters and binary demographics for the open clusters (OCs) NGC 2168 (M35), NGC 7789, NGC 6819, NGC 2682 (M67), NGC 188, and NGC 6791. These clusters span in age from ∼200 Myr (NGC 2168) to more than ∼8 Gyr (NGC 6791) and have all been extensively studied in the literature. We use the Bayesian Analysis of Stellar Evolution software suite to derive the age, distance, reddening, metallicity, binary fraction, and binary mass-ratio posterior distributions for each cluster. We perform a careful analysis of our completeness and also compare our results to previous spectroscopic surveys. For our sample of main-sequence stars with masses between 0.6 and 1M_⊙, we find that these OCs have similar binary fractions that are also broadly consistent with the field multiplicity fraction. Within the clusters, the binary fraction increases dramatically toward the cluster centers, likely a result of mass segregation. Furthermore nearly all clusters show evidence of mass segregation within the single and binary populations. The OC binary fraction increases significantly with cluster age in our sample, possibly due to a combination of mass-segregation and cluster-dissolution processes. We also find a hint of an anticorrelation between binary fraction and cluster central density as well as total cluster mass, possibly due to an increasing frequency of higher-energy close stellar encounters that inhibit long-period binary survival and/or formation.

M35 is a young open cluster and home to an extensive binary population. Using Gaia Data Release 3, Pan-STARRS, and Two Micron All Sky Survey photometry with the Bayesian statistical software, BASE-9, we derive precise cluster parameters, identify single and binary cluster members, and extract their masses. We identify 571 binaries down to GaiaG= 20.3 and a lower limit on the binary frequency off_b= 0.41 ± 0.02. We extend the binary demographics by many magnitudes faint-ward of previous (radial-velocity) studies of this cluster and further away from the cluster center (1.°78, roughly 10 core radii). We find the binary stars to be more centrally concentrated than the single stars in the cluster. Furthermore, we find strong evidence for mass segregation within the binary population itself, with progressively more-massive binary samples becoming more and more centrally concentrated. For the single stars, we find weaker evidence for mass segregation; only the most massive single stars (>2.5M_⊙) appear more centrally concentrated. Given the cluster age of ∼200 Myr, and our derived half-mass relaxation time for the cluster of 230 ± 84 Myr, we estimate ∼47% of the binary stars and ∼12% of single stars in the cluster have had time to become dynamically mass segregated. Importantly, when we investigate only stars with mass segregation timescales greater than the cluster age, we still find the binaries to be more centrally concentrated than the singles, suggesting the binaries may have formed with a primordially different spatial distribution from that of the single stars.

We study the effects of general relativity (GR) on the evolution and alignment of circumbinary disks around binaries on all scales. We implement relativistic apsidal precession of the binary into the hydrodynamics codephantom. We find that the effects of GR can suppress the stable polar alignment of a circumbinary disk, depending on how the relativistic binary apsidal precession timescale compares to the disk nodal precession timescale. Studies of circumbinary disk evolution typically ignore the effects of GR, which is an appropriate simplification for low-mass or widely separated binary systems. In this case, polar alignment occurs, provided that the disks initial misalignment is sufficiently large. However, systems with a very short relativistic precession timescale cannot polar align and instead move toward coplanar alignment. In the intermediate regime where the timescales are similar, the outcome depends upon the properties of the disk. Polar alignment is more likely in the wavelike disk regime (where the disk viscosity parameter is less than the aspect ratio,α<H/r), since the disk is in good radial communication. In the viscous disk regime, disk breaking is more likely. Multiple rings can destructively interact with one another, resulting in short disk lifetimes and the disk moving toward coplanar alignment. Around main-sequence star or stellar mass black hole binaries, polar alignment may be suppressed far from the binary, but in general, the inner parts of the disk can align to polar. Polar alignment may be completely suppressed for disks around supermassive black holes for close binary separations.

With hydrodynamical simulations we examine the evolution of a highly misaligned circumbinary disc around a black hole binary including the effects of general relativity. We show that a disc mass of just a few per cent of the binary mass can significantly increase the binary eccentricity through von-Zeipel–Kozai–Lidov (ZKL) like oscillations provided that the disc lifetime is longer than the ZKL oscillation time-scale. The disc begins as a relatively narrow ring of material far from the binary and spreads radially. When the binary becomes highly eccentric, disc breaking forms an inner disc ring that quickly aligns to polar. The polar ring drives fast retrograde apsidal precession of the binary that weakens the ZKL effect. This allows the binary eccentricity to remain at a high level and may significantly shorten the black hole merger time. The mechanism requires the initial disc inclination relative to the binary to be closer to retrograde than to prograde.

We study the formation of the TRAPPIST-1 (T1) planets starting shortly after Moon-sized bodies form just exterior to the ice line. Our model includes mass growth from pebble accretion and mergers, fragmentation, type-I migration, and eccentricity and inclination dampening from gas drag. We follow the composition evolution of the planets fed by a dust condensation code that tracks how various dust species condense out of the disc as it cools. We use the final planet compositions to calculate the resulting radii of the planets using a new planet interior structure code and explore various interior structure models. Our model reproduces the broader architecture of the T1 system and constrains the initial water mass fraction of the early embryos and the final relative abundances of the major refractory elements. We find that the inner two planets likely experienced giant impacts and fragments from collisions between planetary embryos often seed the small planets that subsequently grow through pebble accretion. Using our composition constraints, we find solutions for a two-layer model, a planet comprised of only a core and mantle, that match observed bulk densities for the two inner planets b and c. This, along with the high number of giant impacts the inner planets experienced, is consistent with recent observations that these planets are likely desiccated. However, two-layer models seem unlikely for most of the remaining outer planets, which suggests that these planets have a primordial hydrosphere. Our composition constraints also indicate that no planets are consistent with a core-free interior structure.

Mutually misaligned circumbinary planets may form in a warped or broken gas disk or from later planet–planet interactions. With numerical simulations and analytic estimates we explore the dynamics of two circumbinary planets with a large mutual inclination. A coplanar inner planet causes prograde apsidal precession of the binary and the stationary inclination for the outer planet is higher for larger outer planet orbital radius. In this case a coplanar outer planet always remains coplanar. On the other hand, a polar inner planet causes retrograde apsidal precession of the binary orbit and the stationary inclination is smaller for larger outer planet orbital radius. For a range of outer planet semimajor axes, an initially coplanar orbit is librating meaning that the outer planet undergoes large tilt oscillations. Circumbinary planets that are highly inclined to the binary are difficult to detect—it is unlikely for a planet to have an inclination below the transit detection limit in the presence of a polar inner planet. These results suggest that there could be a population of circumbinary planets that are undergoing large tilt oscillations.

While giant planet occurrence rates increase with stellar mass, occurrence rates of close-in super-Earths decrease. This is in contradiction to the expectation that the total mass of the planets in a system scale with the protoplanetary disc mass and hence the stellar mass. Since the snow line plays an important role in the planet formation process, we examine differences in the temperature structure of protoplanetary gas discs around stars of different mass. Protoplanetary discs likely contain a dead zone at the mid-plane that is sufficiently cold and dense for the magneto-rotational instability to be suppressed. As material builds up, the outer parts of the dead zone may be heated by self-gravity. The temperature in the disc can be below the snow line temperature far from the star and in the inner parts of a dead zone. The inner icy region has a larger radial extent around smaller mass stars. The increased mass of solid icy material may allow for the in situ formation of larger and more numerous planets close to a low-mass star. Super-Earths that form in the inner icy region may have a composition that includes a significant fraction of volatiles.

We present a fragmentation module and a composition tracking code for the N-body code rebound. Our fragmentation code utilizes previous semi-analytical models and follows an implementation method similar to fragmentation for the N-body code mercury. In our N-body simulations with fragmentation, we decrease the collision and planet formation time-scales by inflating the particle radii by an expansion factor f and experiment with various values of f to understand how expansion factors affect the collision history and final planetary system. As the expansion factor increases, so do the rate of mergers which produces planetary systems with more planets and planets at larger orbits. Additionally, we present a composition tracking code which follows the compositional change of homogeneous bodies as a function of mass exchange and use it to study how fragmentation and the use of an expansion factor affects volatile delivery to the inner terrestrial disc. We find that fragmentation enhances radial mixing relative to perfect merging and that on average, as f increases so does the average water mass fraction of the planets. Radial mixing decreases with increasing f as collisions happen early on, before the bodies have time to grow to excited orbits and move away from their original location.